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ABSTRACT 

We take the line of sight velocity dispersions as functions of radius for 8 Milky Way 
dwarf spheroidal galaxies and use Jeans analysis to calculate the mass-to-light ratios 
(M/L) in Modified Newtonian Dynamics (MOND). Using the latest structural parame- 
ters, distances and variable velocity anisotropy, we find 6/8 dwarfs have sensible M/L 
using only the stellar populations. Sextans and Draco, however, have M/L=9.2^3 q 

and 43.91^9 3 respectively, which poses a problem. Apart from the need for Sextans' 
integrated magnitude to be reviewed, we propose tidal effects intrinsic to MOND, 
testable with numerical simulations, but fully orbit dependant, which are disrupting 
Draco. The creation of the Magellanic Stream is also re-addressed in MOND, the 
scenario being the stream is ram pressure stripped from the SMC as it crosses the 
LMC. 

Key words: gravitation - dark matter - galaxies: clusters 



1 INTRODUCTION 

In addition to the obvious Magellanic Clouds and the 
tidally disrupting Sagittarius dwarf (Ibata, Gilmore & Irwin 
1994), the Milky Way (MW) has 8 relatively nearby dwarf 
spheroidal satellite galaxies for which the radial velocities 
of many hundreds of member stars have been meticulously 
measured at various projected radii allowing the calculation 
of the line of sight (los) velocity dispersion as a function 
of projected radius. If we assume the satellites are spherical 
and in equilibrium, we can use specific mass models to check 
the consistency of the luminous matter in replicating the los 
dispersion profiles through the Jeans equation. 

It so happens that these small satellites appear to be 
some of the most heavily cold dark matter (CDM) dom- 
inated systems in the universe and so including only the 
luminous matter in Jeans analysis underestimates the los 
velocity dispersions considerably. It is because of the well 
studied, small spacial scales and large dark matter content 
that Gilmore et al. (2007ab) have used the relatively low 
dark matter density in the cores of dwarf spheroidals to ar- 
gue for an upper limit on the mass of the DM particles from 
supersymmetric theories. 

There has been speculation that Ursa Minor and Draco 
are formed from the tidal debris of the Magellanic Stream 
(Lynden-Bell 1983; Metz & Kroupa 2007) and that Leo I, 
Leo II and Sculptor have similarly formed from debris cast 
by the tidal break-up of a previously larger Fornax satellite. 
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More hkely, the dwarfs were part of two small groups of 
galaxies that have entered the Milky Way at late times. 

Dwarf galaxies often form from the debris of tidally dis- 
rupted progenitor galaxies. In this case, very little CDM ac- 
companies the stars and gas for the simple reason that the 
tidal disruption must have torn the gas and the CDM out 
into streams from which the gas has condensed and formed 
stars and only a limited amount of CDM can have condensed 
with it. This is contrary to the standard picture where gas 
cools and falls onto an already formed CDM halo. 

There has been a recent study by Bournaud et al. 
(2007) of tidal dwarf galaxies formed from the debris of 
NGC5291. They found a persisting mass discrepancy by 
measuring the visible galactic masses and comparing with 
the mass implied by the rotation curves. They proposed 
the existence of hitherto unheard of disk DM, in the pro- 
genitor galaxies, in the form of molecular gas (as opposed 
to the cold DM of massive particle form). Two simulta- 
neous studies (Milgrom 2007 and Gentile et al. 2007) in 
the framework of Milgrom (1983) 's Modified Newtonian Dy- 
namics (MOND: for reviews see Sanders & McGaugh 2002; 
Bekenstein 2006; Milgrom 2008) found the discrepancy was 
exactly that expected in MOND because the tidal dwarfs 
are very diffuse and low surface brightness (McGaugh & 
de Blok 1998) meaning they have low internal accelerations 
(much less than the transition acceleration of the theory 
ao = 1.2 X 10""ms"^ = 3.6( kms~^)^pc'\ In addition, 
the external accelerations from the large progenitors is also 
weak (discussed in depth in both papers) allowing a rela- 
tively large boost to the internal gravity enabling perfect 
agreement with the observed rotation curve using just the 
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luminous matter (no need for CDM or molecular disk DM). 
This is the most remarkable evidence for MOND since the 
baryonic TuUy-Fisher relation was matched in great detail 
(McGaugh et al. 2000; McGaugh 2005). 

Previous attempts to check the consistency of the Milky 
Way's dwarfs with MOND have been limited to single num- 
ber central velocity dispersions (although see Lokas et al. 
2001 which sadly does not include the external field ef- 
fect). This gave good estimates of the mass-to-light ratios 
(M/L) for all the dwarfs studied here by considering whether 
each dwarf was dominated gravitationally by the external 
field of the Milky Way (MW) or by its own internal grav- 
ity (see Milgrom 1995, Brada & Milgrom 2000). If it was 
dominated by the external gravity, g^x then the Newtonian 
M/L was used after renormalisation of Newton's constant, 
Geff/G = aa/gex, reducing the M/L in MOND by the same 
factor. If the external gravity was weak w.r.t. the internal 
gravity then the dwarf mass was calculated from the re- 

lation M — xsf^ from which the M/L can be deduced. 
Nevertheless, proper Jeans modelling with the luminosity 
profile of the dwarf along with the near exact calculation of 
the dwarf's internal gravity at all radii by including analy- 
sis of the external field effect (see Milgrom 1983, Angus & 
McGaugh 2008ab and Wu et al. 2007) and variation of the 
velocity anisotropy profile (Angus et al. 2008) in order to 
match the observed los dispersions at all radii is far more 
satisfactory. This also gives tighter limits on the M/L. 



2 MODELING THE LOS VELOCITY 
DISPERSIONS 

The key data for the analysis are the luminosity of the dwarf, 
its structural profile and the Galactocentric distance. The 
total luminosities, L„ have been measured only in the V- 
band and some, especially that of Sextans, are very old and 
we suggest that new K-band images would be more appro- 
priate measures of the stellar mass (Driver et al. 2007). 

Next we use the King models (concentration c and tidal 
radius Rt) for the enclosed luminosity profiles, Lxir), of the 
dwarfs (mostly taken from Irwin & Hatzidimitriou 1995) 
that allows us to numerically calculate the luminosity den- 
sity as a function of radius. Finally, we need the current 
Galactocentric radius, Rmw in order to compute the exter- 
nal gravity from the MW via the flat rotation curve (see 
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amplitude of the rotation curve is from the SDSS data com- 
piled by Xue et al. (2008) . From this we solve for the internal 
gravity 



g{r)^i{\g(r) + gex\/ao) = G{M / L)LK{r)r- 



(1) 



where we use the function employed by Famaey & Binney 
(2005) in their fit to the MW's terminal velocity curve 



^i(x) = x/{l + x), 



(2) 



and we take \g{r)+ g^x \ = a/ g('i')'^ + gix which only matters 
in the small transition region. Then we numerically differ- 
entiate the luminosity profile from the King model, Ljf (r), 
to match with the simple function 



l(r) oc 1 + — 



which is a luminosity density profile and helps us to analyt- 
ically define a[r) = d\nl[r) / d\nr — —Oo for solving 
Jeans equation for the radial velocity dispersions, ar, 



d_ 

dr 



ar{r) + 



Q(r)+2/3(r) , 



<^r{r) = -g{r) 



(4) 



and we can cast this into the los with the use of the lumi- 
nosity density profile, l{r) (see Binney & Mamon 1982) as 
per Angus et al. (2008). Taking an anisotropy, /3 of the form 



f3{r)=Po-—. 



(5) 



This and all other fit parameters are given in table 1. 

Here we have three free parameters available for match- 
ing the observed and model velocity dispersions. They are 
the M/L and the two associated with /3. M/L adjusts the 
normalisation, whereas /9 varies the shape of the los pro- 
file. Note that although there are two parameters for the 
anisotropy, the profile required is virtually the same for 6/8 
dwarfs (see fig 2). In any case, Koch et al. (2007) make it 
clear from their analysis of Leo II that there is compelling 
evidence for velocity anisotropy in dwarfs. 

By fitting the los velocity dispersions of Mateo, Ol- 
szewski & Walker (2007) for Leo I, Muiioz et al. (2005) for 
Ursa Minor and Walker et al. (2007) for the remaining 6, 
with the model MOND dispersions, we were able to best fit 
the anisotropy profile of all 8 dwarfs and thereby compute 
the M/L. The fits are not hugely sensitive to the anisotropy 
profile. For instance, in table 1 the column of also gives 
the 1-(T errors on that parameter (found by simply hold- 
ing M/L and /3o constant and varying . Most dwarfs have 
very little sensitivity to the rate at which the anisotropy de- 
creases with radius, with the only exception being Draco. 
Many of the dwarfs (Leo I, Leo II, UMi and Sextans) are 
quite comfortable with constant anistropies, however, none 
of the dwarfs are consistent with an anisotropy that becomes 
highly biased towards tangential orbits in small increases of 
radius. It is difficult to be quantitative about how we could 
rule out MOND with any of these velocity anisotropies, since 
we use a very simple model for this and anisotropy could 
vary in a non-linear way, or there could be distinct stel- 
lar populations with completely different anisotropies. How- 
ever, if instead of becoming biased towards tangential orbits 
at large radius, they were in reality biased towards highly 
radial orbits, there would be a serious problem because no 
M/L without a dark halo could fit the shape of the profiles. 



where the An example of this is given for the Sculptor dwarf in Fig 1(f) 



(3) 



where the red line shows the fit with a constant /3(r) = 0.8. 
No scaling with M/L can produce a good fit and it shows 
how precise knowledge of the velocity anisotropy would be 
extremely constraining to MOND, whereas any CDM halo 
can never be constrained to such detail and thus can never 
be falsified. 



3 RESULTS 

The fits to the los velocity dispersion profiles of all 8 dwarfs 
are shown in Fig 1. The solid line is the best fit and the 
M/L is given in table 1 along with all other relevant data 
for each dwarf. The dashed lines are the 1-cr errors in M/L. 
The reduced /{n — Ud.f ~ 1), where n is the number of 
data points and Ud.f is the number of degrees of freedom in 
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Name 




c 


Rt 




Lv 


Vr 


rex 






l3o 


M/L„ 




^red 


Carina 


kpc 
101 ± 5 


0.51±0.08 


PC 

581 ± 86 


-9.3±0.2 


4.4±1.1 


km s~ "'^ 
8 


PC 

2000 


24 


PC 


0.6 




80 


0.72 


Draco 


93 ± 6 


0.72±0.05 


1225 ± 80 


-9.0±0.3 


3.3±1.1 


-294 


1500 


14 




0.93 


--^?r3 


270 


0.65 


Leo I 


257 ± 8 


0.39±0.07 


1002 ± 50 


-12.4±0.2 


76±7 


178 


3700 


22 




0.1 


„ -+0.65 
-0.3 


10 


0.47 




95.5 ± 3 


n.98±n.i4 


3445 ± 1141 


-9.7±0.2 


6.3±1.4 


75 


1200 


6 




0.0 




102 


1.2 


Fornax 


138 ± 8 


n.72±n.05 


2078 ± 177 


-13, 2 ±0,2 


158±16 


-36 


3200 


14 


„„+1500 


0.3 


1 4 + 0.45 
■ -0.35 


11 


1.4 


Sculptor 


87 ± 4 


i.i2±n.i2 


1329 ± 107 


-11,1 ±0,2 


28±4 


95 


1200 


14 




0,5 




36 


0.7 


Leo II 


233 ± 15 


0.48±0.1 


554 ± 68 


-9.9±0.3 


7.6±3 


22 


2500 


30 




0.56 


1.85+1° 


57 


1.1 


UMi 


76 ± 4 


0.64±0.05 


1977 ± 104 


-10.3±0.4 


11.0±4.8 


-87 


6000 


25 




0.73 




440 


0.86 



Table 1. All distance and magnitude data from Mateo (1998), Galactocentric radial velocities from Wilkinson & Evans (1999) and 
structural parameters from Irwin & Hatzidimitriou (1995) except Draco's Vr which is from Hargrcaves ct al. (1996). Rmw actually 
Rq, but since satellites have Galactic latitudes such that 90° < I < 270°, this will only underestimate Rmw- Recalculated distance 
(Bellaazini et al. 2002 for both) and structural values of are used for Draco (Odenkirchen et al. 2001) and Ursa Minor (Palma et al. 2003) 
which lead to new absolute magnitudes. Sextans has an updated distance from Lee et al. (2003) which gives a 20% increased luminosity. 
Leo Fs luminosity and structural parameters come from Sohn et al. (2007) which makes use of the new distance modulus from Bellazzini 
et al. (2004). Luminosity for Leo II corrected for new distance from Bellazzini et al. (2005), as is Sculptor's from Piatek et al. (2006). The 
solar absolute magnitude is Mt, = 4.8. Wo parameters for the King model are 1.97, 3.14, 1.51, 4.59, 3.14, 5.29, 1.85, 2.7. The Newtonian 
M/L come from dividing the dark mass at Tmax from Walker et al. (2007) by the V-band luminosity (except UMi which comes from 
Kleyna et al. 2002). n is the number of data points and n^./ is the number of degrees of freedom in the model (always 3). 



the model (always 3), is very reasonable for each dwarf, even 
though the data are inherently noisy. A good example of this 
is the Fornax dwarf in which there are a larger number of 
data points (15) at large radii (from 500-lOOOpc) for which 
the individual errors are relatively small compared to the 
spread in the data. No straight (or even with curvature in 
one direction) line through these data points could obtain 
a low reduced and the line of best fit to these data is a 
very good one. 

The anisotropy profiles for each dwarf are plotted in 
Fig 2 against R/Rc and the linotype is given in tabic 1. 
It is interesting that, apart from Ursa Minor and Draco, 
the other 6 dwarfs have very similar anisotropies going from 
isotropic or slightly radially biased orbits at their centres to 
more tangcntially biased orbits at the edges. 

To be comfortably explained in MOND, these dwarfs 
should have M/L of around 3 solar units. It is good to see 
that all four dwarfs associated with the Fornax- Leo-Sculptor 
stream have normal M/L. This is absolutely necessary be- 
cause Leo I and Fornax are very massive compared to the 
rest which means the external gravity of the MW should not 
disturb them. Furthermore, the two Leo dwarfs are beyond 
200kpc meaning even if they had a near pericentre to the 
MW, they should have returned to equilibrium by now. 

The Carina and Ursa Minor dwarfs have borderline rea- 
sonable M/L (5.6^2.9 and S.Stg J respectively). When the 
l-fT error in M/L is considered along with deviations from 
spherical symmetry, coupled with other uncertainties in lu- 
minosity, distance, structure and possible tidal effects, they 
may well be safely within the bounds of normal M/L. 

Although initially thought to have a large mass discrep- 
ancy, reanalysis of Ursa Minor's structural parameters and 
distance by Bellazzini et al. (2002) showed it to be more 
distant and extended, therefore it has a significantly larger 
luminosity which has reduced the M/L. 



3.1 Problem cases 

There can be little doubt that Draco and now Sextans are 
both concerns for MOND. Their M/L's are uncomfortably 
high and we seek to find an explanation for this given that 



the other 6 are more or less compatible. First of all, it is 
troubling that these two have the noisiest velocity disper- 
sions especially compared to the almost constant velocity 
dispersions of Fornax, Sculptor, Leo I and even Carina. The 
dispersions seem to have regions where they are abnormally 
high and low, in addition to the presense of substructure in 
Sextans (Kleyna et al. 2004). Nevertheless, even if we only 
paid attention to the low values (i.e. with the lower dashed 
lines in Fig 1, the M/L are still unsatisfactory. Draco's lu- 
minosity and distance have recently been reviewed by Bel- 
lazzini ct al. (2002), but Sextans still has a very poorly con- 
strained structure and it would be of great benefit to have 
its King model parameters and integrated magnitude re- 
analysed. Still, it is worth bearing in mind that although 
the M/L are high, they are also very poorly constrained and 
they are only ~2-cr from acceptable values. 



3.2 The double impact of Galactocentric distance 
in MOND 

For a satellite like any one of the dwarfs studied here, the 
only infiuence the Newtonian gravity of the MW exerts is to 
carve out its orbital path. Any tidal forces are negligible and 
there is absolutely no difference to the internal gravity. In 
stark contrast, the MW's MOND gravity has a huge impact 
on the internal gravity of the low luminosity, nearby dwarfs: 
Carina, Sextans, Draco and Ursa Minor. As mentioned 
above, the external field effect renormalises the effective 
Newtonian Gravitational constant Geff = G/^{ ^ ^™ — ) (as 
long as the external gravity is higher than the internafgrav- 
ity). As Rmw increases, Ge// follows suit reducing all in- 
ferred M/L by the same factor. The key thing here is that 
if new observations increase the distance to a dwarf by 11% 
as Lee ct al. (2003) has done for Sextans (and Bellazzini ct 
al. 2002, 2004, 2005 have done for other dwarfs) not only 
does the luminosity increase by 21%, which obviously drops 
any M/L by 17%, but the reduction in the external field 
effect allows the M/L to drop a further 10% (in that partic- 
ular case). So, the increased distance of less than lOkpc for 
Sextans dropped the M/L from 15 to almost 9. 
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3.3 Sextans: poorly constrained light profile 

Sextans has a M/L=9.2l3 g which is impossible from popu- 
lation synthesis point of view. However, this M/L assumes 
we know the luminosity of Sextans well. In fact, Sextans is 
the lowest surface brightness dwarf we study here and it's 
light profile has not been examined in detail since before 
Mateo (1998). All of the other dwarfs that have had their 
luminosity re-evaluated have had it pushed up far beyond 
the original 1-a error: UMi, Draco, Carina, Leo I and Leo 
II. In fact, the first four of those five mentioned saw their 
luminosities double. When Sextans has its luminosity re- 
integrated with new photometry, we expect it to rise by at 
least 50%, but not necessarily by a factor of 3 which would 
put it in the acceptable zone. This will probably leave it 
with a similar M/L as Carina and UMi, requiring a modest 
amount of tidal heating (see next section) to be consistent 
with the los velocity dispersions. 

At the moment however, it should be noted that al- 
though the M/L of Carina and UMi are higher than 3, un- 
certainties in their los velocity dispersions as well as their 
luminosity and distance mean they are in agreement within 
the observational errors. Still, the fact that they are all co- 
herently, higher than preferred compells us to consider why 
in the next section. 



3.4 Impulsive tidal heating 

The most enticing resolution of the three moderately high 
M/L dwarfs are that Sextans, Carina and UMi have recently 
had a close pericentre of the MW and have been slightly 
tidally heated. The crossing times {2Yi/aios) at the last mea- 
sured los dispersion radius are 120 (Carina), 300 (Sextans), 
150 (Umi) and 110 (Draco) Myr, but recall that relaxation 
times in MOND are higher than in Newtonian gravity as 
calculated from galaxy mergers by Nipoti et al. (2007) and 
Tiret & Combes (2007). Using a simple back stepping orbit 
integrator taking the current Galactocentric radii and ra- 
dial velocities, with gravity calculated from a flat rotation 
curve with circular speed of 170 kms~^ (see Xue et al. 2008), 
we found it has been roughly 700 (Carina), 500 (Sextans) 
and 850 (UMi) Myr since the satellites had their last peri- 
centre. Draco's apocentre is over 400kpc, which means it is 
probably on its first passage through the MW since this or- 
bit integration neglects the external field of M31 and large 
scale structure on the MW and despite this it is 5600Myr 
since Draco left. 

There is something key to remember here in which 
MOND differs starkly from Newtonian gravity. A low surface 
brightness dwarf galaxy has dynamical stability imparted 
by the MOND enhanced gravity, much like the DM halo 
provides for Newtonian gravity. If we make the crude but 
fair assumption, for clarity's sake, that the external field 
of the MW dominates the internal accelerations of the 3 
satellites, then the MOND effect is to renormalise Newton's 
constant upwards which enhances the internal gravity. So 
Ge// = G/ ^{gex/ao) ~ Gao/g^x- 

As the satellites approach the MW this stability is lost 
due to the rising external field which bleeds the internal 
gravity of the dwarf. This of course doesn't happen in New- 
tonian gravity and the dwarfs can withstand tidal forces 
more easily (the DM halo stays intact). This makes it far 



easier to perturb a MOND dwarf. The radii below which 
tidal effects begin to affect the dwarfs ~40kpc as opposed 
to <10kpc for Newtonian gravity with the observed DM ha- 
los (see Brada & Milgrom 2000b). 

However, currently, none of the dwarfs studied here 
are strictly within the range of radii with measured disper- 
sions, in the non-adiabatic/impulsive regime. However, they 
may have just left this region. We re-emphasise their point 
that tidal effects begin to impose themselves at much larger 
Galactocentric radii in MOND. We slightly repackage this 
idea by looking at the rate of change of the external field 
strength during one crossing time at the core radius of each 
dwarf i.e. ± (^) ~ £2L _ plotted 
against M/L for each dwarf in Fig[T}. 

From Fig 3 one can see that low Galactocentric ra- 
dius, luminosity and surface brightness dwarfs have higher 
M/L. This is unsurprising because these properties make the 
dwarf more susceptable to tidal effects. High rate of change 
of external field is where Draco stands out with its enor- 
mous M/L. Although the dwarfs are not in the non-adiabatic 
regime, the magnitude of the external field will have been 
varying rapidly in the past. It is unknown at present how 
serious an effect this is, but it needs to be studied with re- 
alistic, high resolution MOND simulations (see Nipoti et al. 
2007; Tiret & Combes 2007) for dwarf galaxies in the exter- 
nal field of the MW, but also for spirals passing through the 
centres of clusters. This rapidly varying external field we pre- 
dict is destroying the integrity of spirals in clusters (see Ru- 
bin et al. 1999) and also imposes a morphology-position rela- 
tion in clusters more easily than galaxy harassment (Moore 
1999). 

Whether Sextans or any of the others have recently had 
a close pericentre passage will only be known after accurate 
proper motions are published. The proper motions of Ursa 
Minor are known well enough only to say the pericentre 
lies between lOkpc and 76kpc (Piatek et al. 2005) at 1-cr. 
Brada & Milgrom (2000) showed that before tidal effects 
are felt there is a region where non-adiabatic (impulsive) 
effects due to the quickly changing external field heat up 
the dwarf without breaching it. There may be insufficient 
time for the satellite to regain equilibrium after leaving this 
zone, especially since they are still external field dominated. 

Of course, when accurate proper motions become avail- 
able, if the opposite is true and they are currently at peri- 
centre and on a reasonably circular orbit, they become quite 
a serious problem for MOND. Finally, as mentioned in i]3.1l 
Sextans should be reanalysed in more detail to tighten the 
errors on its light profile and magnitude. 



3.5 Draco: tidal disruption 

If indeed Sextans has been heated by a close pericentre, the 
biggest issue we are left with is Draco with M/L=43.9l^g'3. 
Draco is now moving towards the MW at ~ 294kms~^ 
which means it is most likely on its first passage of the MW. 
Segall et al. (2007) showed that Draco has a symmetric, un- 
perturbed density profile with no evidence for tidal tails - 
"Draco has led a quiet existence". This is in line with the 
fact that it has been on a direct path towards the MW since 
its creation. However, if its tangential velocity is negligible, 
its tidal tails may line up along the los. 
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Since Draco cannot have been tidally heated by a 
near pericentre passage due to the timescales involved, this 
leaves us with only one conclusion i.e. that Draco is in 
the throes of tidal disruption, which is perfectly possible 
in MOND, but unimaginable in Newtonian gravity. At its 
apocentre, Draco would have experienced very little exter- 
nal field and thus its low surface density would be enough 
to maintain its dynamical integrity due to the MOND en- 
hanced gravity. This would make its internal gravity at 
R^lOOOpc ~ ^3.6 X 4.4 x IQ-^ x 3 x 3.3 x lOVlOOO^ ~ 
ao/30 (this assumes no external field). However, currently 
it is at 93kpc and thus the internal gravity is merely 
4.4x10"^ X ao^^^§— ~ ao/70, reduced because of the grow- 
ing external field. This loss in binding energy might have 
recently thrown Draco into the early stages of tidal disrup- 
tion. Munoz et al. (2007) have showed with high resolution 
N-body simulations that even modest tidal disruption near 
the edges of dwarfs can have a dramatic impact on the ob- 
served dynamics. This again must be confirmed by numer- 
ical simulations of Draco on a plunging (radial) orbit from 
a large radius in MOND with a reasonable M/L and veloc- 
ity dispersion found by solving the Jeans equation. Either 
Draco will begin to disrupt at ~ lOOfcpc and its los velocity 
dispersions will appear artificially high, or it will have dis- 
persions consistent with the Jeans equation at all radii down 
to at least 50kpc. 



3.6 Dwarf fly-bys 

Is it possible that the dwarfs are co-responsible for their 
own large dispersions? Draco and Ursa Minor are sep- 
arated by approximately 20kpc, so the largest exter- 
nal gravity that UMi could inflict on Draco is a mere 
Gao/gex{M/L^)Lsatd~^ ~ 4.4 x 10"^ x 10 x 3 x I.l x 10*' x 
(20000)-^ ~ flo/lOOOO, and since they are moving towards 
each other now, they were never any closer in the past, so 
there is no chance of this being important. If say the Fornax 
satellite conspired to collide with Sextans it would have to 
pass within d ~ ,/G (M/L^)Lsat ~ V4.4 x 10-^ ~ 
4500j3C. That would balance the internal and external ac- 
celerations because at Ikpc the internal gravity of Sextans 
is ao/20. This is at least possible given that they are cur- 
rently moving away from each other. Accurate proper mo- 
tions would resolve this. 



3.7 Tidal tails and ellipticity 

It was noted by Klimentowski et al. (2007) and Lokas et 
al. (2008), tidal tails along the line of sight can have an ad- 
verse effect on the los velocity dispersion. They showed using 
high resolution N-body simulations as comparisons that the 
velocity dispersion is boosted by the interlopers, giving the 
appearance of a high M/L. This may indeed reduce the M/L 
for several of the dwarf spheroidals studied here, but it re- 
mains to be seen whether the effect is signiflcant enough to 
reduce Draco's M/L by an order of magnitude. 

Linked to this is the idea that deviations from spherical 
symmetry could produce the high dispersions. As above, if 
the major axis was pointing along the los, then this would 
have the same effect as tidal tails. So significant help could 



be granted by tidal streams or deviations from spherical 
symmetry. 

3.8 Relativistic MOND 

The field of relativistic MOND kicked off by Bekenstein 
(2004) is growing fast and the best bet at the moment is the 
Einstein Aether theory (see Zlosnik, Ferriera & Starkman 
2007), however, one should see Milgrom (2008) for a review 
of possible fundamental origins for MOND. Currently, it is 
unknown exactly how systems, like dwarfs, with rapidly time 
varying external gravity might fare in these theories, but it 
may be the case that the vector-field boosts the internal 
gravity to a greater extent than simple MOND. 

3.9 Dust 

The easiest solution to the problem would of course be if 
there was simply a large dust sheet obscuring most of the 
light of the three galaxies. In a recent paper by Driver et 
al. (2007) they postulated that more than 90% of the light 
in galaxy bulges (dwarfs are many orders of magnitude re- 
moved bulges in most properties) can be obscured in the 
V-band by dust. This effect actually is still not enough and 
the remaining light would have to be attenuated by dust 
in the Milky Way. Draco and Ursa Minor are on the side 
of the Milky way opposite the bulge from the solar posi- 
tion, so it is not even as if we are looking through high 
stellar, dust and gas density towards them. Finally, since 
deep Color-Magnitude-Diagrams exist for all these galaxies, 
if they were affected by significantly larger extinction than 
presently believed it would have been noted long ago. 

3.10 Trends with key variables 

In Fig|3]we plot the M/L against four variables which de- 
scribe the dwarfs in order to elucidate trends, although we 
are clearly dealing with large errors and small number statis- 
tics. In the top left panel, we plot M/L against Galactocen- 
tric distance. There is an obvious division between Fornax 
and the two Leo dwarfs beyond 130kpc and the other 5 less 
than lOOkpc from the centre of the MW. In addition, there 
is some evidence for a trend shown in the top right panel for 
lower surface brightness dwarfs to have higher M/L as well 
as those with an external field that is varying more rapidly. 

However, none of those three independent variables cut 
a tight correlation with M/L because none of them are the 
cause of the high M/L, they are merely strong indicators. 
The smaller the Galactocentric radius, the higher the exter- 
nal field and the more energy there is available to perturb 
the dwarf. The lower the luminosity, the more susceptible the 
dwarf is to being perturbed by the external field. The radial 
velocity (bottom left) is a poor indicator because any per- 
turbing of the dwarfs was done near pericentre. This again 
shows why Galactocentric radius is just an indicator because 
it simply means the dwarf was more likely to have been at 
smaller radii but for all we know it may currently be at peri- 
centre (Bellazzini et al. 1996). It is not surprising that Draco 
with its large M/L has a low Galactocentric radius, surface 
brightness and large velocity towards the MW as it may well 
be tidally disrupting. 
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Figure 1. Velocity dispersion profiles, cr;oa(R), for all 8 dwarfs. Best fit lines have solid line type and the dashed curves correspond to 
the 1-sigma errors on M/L. 



4 THE MAGELLANIC STREAM 

The new proper motions of the Magellanic Clouds (Kalli- 
vayalil et al. 2006ab) show conclusively that the Large Mag- 
ellanic Cloud (LMC) and Small Magellanic Cloud (SMC) 
binary are currently at pericentre (Besla et al. 2007, here- 
after B07) in their orbit of the MW. The creation of the 
Magellanic stream, a highly coUimated gas cloud arcing 180° 
across the sky and emanating from the cloud's position. Sev- 
eral puzzling features plague the stream such as the fact the 



stream does not trace the orbit of the LMC in any favourable 
CDM (B07, Fig. 9) or MOND orbit (See X.F. Wu et al. in 
prep), there are very few stars in the cloud (B07) and at a 
pericentre of > 50kpc the stream cannot be tidally induced. 
Eyeballing Fig. 18 of B07 one can see the motion of the SMC 
is significantly influenced by the LMC allowing its orbit to 
trace the stream. As to how the gas was stripped from the 
SMC, we need only consider the ram pressure of the gas 
in the LMC. As mentioned above, on close approach to the 
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Figure 3. Shows M/L against several variables for each dSph. (a) Galactocentric radius, which implies external field strength (b) V 
band surface density (c) Galactocentric radial velocity (d) time derivative of the external field in units of the dSph's dynamical time at 
the core radius. Note the three dwarfs with large time derivative of the external field are Draco, Ursa Minor and Sextans. 



LMC, the SMC's internal gravity is damped by both by the 
MW and more importantly, by the LMC. 

The criteria for ram pressure stripping (Roediger & 
Bruggen 2007) is for the gravitational restoring force per 
unit area, 



fgraviR) = T.g ,S MC (R) 9 z{R) , 

to be less than the ram pressure 

Pram = pg,LMcVrel- 



(6) 



(7) 

where Eg^s is the surface density of gas in the SMC, Qz is 
the SMC's gravity in the direction of the collision, pg is the 
density of gas in the LMC and Vrei is the relative velocity 
between the LMC and SMC at collision. 

The relative velocity between the pair at the moment 



is around Vrei = lOOkms"^. The total gas mass of the 

LMC is Mg.LMC = 5 X 10^ M©, so if it has dimensions 
(R,z)=(8.5kpc, ±0.2kpc) then the average density of gas 

is Pg,LMC = ^(8.5x10^^x400 ~ 0.02Mopc"^. The mass of 
gas in the SMC is Mg^sMC = 4.2 x IQ^Mq (Stanimirovic, 
Staveley-Smith and Jones 2004)so for a typical radius of 
2.5kpc, the surface density of gas is T,g^sMC = 22Mopc~^. 

The key input to this equality and the only variable 
that changes from MOND to Newtonian dynamics is the 
gravity of the SMC in the direction of the collision. The 
total baryonic mass of the SMC is around Msmc ~ 8.5 x 
10** M© (Stanimirovic et al. 2004). With the LMC so close 
to it, the SMC's internal gravity is simply Newtonian with 
a renormalised gravitational constant, Ge// = Gao/gLMC- 
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Figure 2. Anisotropy profiles, /3(r), for all 8 dwarfs. See table 1 
for the corresponding linetypes. 

At 3.5kpc, the rotation curve of the LMC is fiat (Piatek, 
Pryor and Olszewski 2007) at 120kms~^. At this radius, 
Qlmc ~ 1.2ao. This means that Geff < 2G, whereas if 
the LMC was not there, and (more or less) compared to 
Newtonian gravity, Geff ~ 7G. 

With these special case parameters, the gravity must be 
(?z ~ 2.5ao to retain the gas, whereas the gravity is less than 
this; gz,M(R = 3.5kpc) ~ ao/6 and the Newtonian value is 
3-4 times larger. This shows that ram pressure stripping in 
MOND is easier which is completely contrary to the naive 
assumption ignoring the external field effect. This can easily 
be the mechanism which removes the gas from the SMC and 
creates the Magellanic Stream. This has further implications 
for the ram pressure stripping of normal spirals in clusters 
and other satellite galaxies of the MW. 



5 CONCLUSION 

Here I have studied the los velocity dispersion profiles of 8 
of the MW's dwarf satellites in MOND. We have shown that 

6 of the dwarfs (Leo I, Leo II, Fornax, Sculptor, Ursa Mi- 
nor and Carina) have acceptable M/L within their errors. 
Despite Carina and Ursa Minor having M/Ls of S.Glj'g a-nd 
5.8l3 g respectively, it is pointless to be too concerned about 
them because they have M/L of 3 within 1-a and uncertain- 
ties in luminosity, ellipticity, structure and equilibrium could 
easily account for this discrepancy. On the other hand, the 
remaining 2 (Draco and Sextans) have M/L completely in- 
consistent with the stellar populations. The possible causes 
of this have been discussed in depth and 3 testable effects 
that cause the discrepancies have been suggested. 

First of all, the photometry of Sextans is extremely old 
(more than 10 years) and so I suggest K-band images of 
the dwarf should be acquired to alleviate any doubts on the 
luminosity of the dwarf and how the light is distributed. 
However, this is likely to leave Sextans with a M/L similar 



to Carina, UMi and Sextans (~ 5). To be certain of the 
effects of non-adiabatic tidal heating of the dwarfs, we need 
accurate proper motions to make sure Sextans is currently 
not at pericentre and that Draco is on a highly elongated 
orbit plunging towards the MW. 

If this is the case, then realistic, high resolution sim- 
ulations (like those of Nipoti et al. 2007 and Tiret & 
Combes 2007) must be run in order to confirm that the non- 
adiabatically changing external field effect can indeed make 
the small boost to the los velocity dispersions in these satel- 
lites over the given orbits and timescales and even tidally 
disrupt Draco. If it turns out that the external field effect 
is unable to disrupt Draco given its orbital history then it 
poses a significant problem for MOND. 

It is quite remarkable that using only the stellar mass 
of the dwarfs, we can successfully match the los velocity 
dispersion profiles with sensible falling anisotropics and get 
excellent reduced for all 8 dwarfs and 6/8 have good 
M/L. As for draco, the unknown orbital history is vital to 
the comprehension of its large M/L and until it is known, we 
cannot confirm the dwarf galaxies of the MW as a success 
or failure of MOND. 
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